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Abstract 

We consider resonant Dirac leptogenesis in a geometry with three five-dimensional throats 
in the flat limit. The baryon asymmetry in the universe is generated by resonant decays of 
heavy Kaluza-Klein scalars that are copies of the standard model Higgs. Discrete exchange 
symmetries between the throats are responsible for establishing two key features of the model. 
First, they ensure a near degeneracy of the scalar masses and thus a resonant decay of the 
scalars. This allows for Dirac leptogenesis at low energies close to the TeV scale. Second, 
the discrete symmetries connect the observed baryon asymmetry with the Yukawa couplings 
of the low-energy theory. As a consequence, we obtain correlations between the low-energy 
leptonic mixing parameters and the Dirac CP phase that can be tested at future neutrino 
oscillation experiments such as neutrino factories. 
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1 Introduction 



One of the major questions in neutrino physics is whether neutrinos are Dirac or Majorana 
particles. Currently, considerable experimental effort is underway [f] to measure neutrino- 
less double beta decay (OvfiP), which would require the neutrinos to be of Majorana-type. 
Neutrino oscillation experiments, however, cannot distinguish between Dirac and Majorana 
neutrinos and as long as Ou/3/3 has not been measured, there will always be the possibility 
that neutrinos are Dirac particles - just like all other fermions in the standard model (SM). 

One advantage of having Majorana neutrinos is that the smallness of the observed light 
neutrino masses ~ 1CT 1 eV could be understood in terms of the seesaw mechanism [2,3] which 
establishes a connection to grand unified theories (GUTs). The type-I seesaw mechanism [2] 
offers, furthermore, the possibility to understand the observed baryon asymmetry in the 
universe (BAU) [4,5] through baryogenesis via leptogenesis [6]. In standard leptogenesis, the 
BAU is produced by the decay of the heavy SM singlet Majorana neutrinos that generate 
small neutrino masses via the seesaw mechanism (for reviews see, e.g., [7] and [8]). If, 
however, the neutrinos are Dirac particles, the original leptogenesis scenario would no longer 
apply since the SM singlet neutrinos would have zero Majorana mass. The BAU can then, 
instead, be generated by Dirac leptogenesis [9]. Several studies have shown that Dirac 
leptogenesis may indeed be responsible for the BAU [10-13]. 

In Dirac leptogenesis, the BAU is generated by the decay of heavy copies of the SM 
Higgs doublet into leptons. Due to the smallness of the Dirac neutrino Yukawa couplings 
the resulting lepton asymmetry can be stored sufficiently long in the singlet neutrino sector 
to allow for successful baryogenesis via sphaleron processes [14]. The original version of Dirac 
leptogenesis, however, raises a couple of questions. First, the scenario does not address the 
origin of the heavy copies of the SM Higgs. Second, Dirac leptogenesis would serve as a 
GUT-scale leptogenesis scenario [15] by preferably taking place at rather large energies near 
~ 10 16 GeV. This may, however, get into conflict with standard inflationary models and the 
gravitino problem [16]. Besides that, the Yukawa couplings responsible for leptogenesis seem 
to be completely unrelated to the Yukawa couplings giving rise to the observed neutrino 
masses. In other words, in the original scenario for Dirac leptogenesis, a measurement of the 
low-energy lepton mass and mixing parameters would have no connection with the BAU. 

In this paper, we consider a model for Dirac leptogenesis that addresses all of these prob- 
lems. The model makes use of discrete symmetries to (i) implement resonant leptogenesis 
at low energies close to the scale of electroweak symmetry breaking (EWSB) [17, 18] (see 
also [19]) and to (ii) relate the observed BAU with the low-energy lepton mixing parameters 
measurable in neutrino oscillation experiments. For this purpose, we work in the flat limit of 
a five-dimensional (5D) background with several "throats" that can emerge from flux com- 
pactification in string theory [20,21]. The field theory on this multi-throat background [22,23] 
allows to identify the heavy scalars necessary for Dirac leptogenesis with Kaluza-Klein (KK) 
excitations and use the field separation on the throats to naturally implement nearly mass 
degenerate scalars that decay resonantly. 

The paper is organized as follows: In Sec. |2j we briefly review the idea of Dirac lepto- 
genesis. Next, in Sec. |3j we present our model for Dirac leptogenesis on a background with 
three throats. In Sec. [4j we discuss the boundary conditions of the 5D scalars and fermions 
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Figure 1: Generation of baryon asymmetry in Dirac leptogenesis [9]. Starting from LH 
and RH lepton asymmetries Lr = —Ll > in the neutrino sector, sphaleron processes (1) 
produce a baryon asymmetry AB = ALl from Ll. For sufficiently small Dirac neutrino 
Yukawa couplings, LR equilibration processes (2) set in later and lead to a final total lepton 
asymmetry L" = L" L + L" R = B" = AB. 



along with the resulting wavefunction profiles and the Yukawa couplings. Then, in Sec. [5j 
we determine the range of Yukawa couplings necessary for successful Dirac leptogenesis and 
discuss the connection of the BAU with the low-energy lepton mixing parameters. In Sec. [6j 
we present our summary and conclusions. Finally, in the appendix, we give further examples 
for correlations of the low-energy mixing parameters. 

2 Brief Review of Dirac Leptogenesis 

Let us start by giving a short review of the Dirac leptogenesis scenario proposed in [9]. We 
assume the SM gauge group Gsm = SU(3)c x SU(2)l x U(1)y- Unless otherwise stated, 
it is, in the following, understood that a left-handed (LH) fermion is in an SU(2) L doublet 
representation while a right-handed (RH) fermion is an SU(2)l singlet. We denote the LH 
lepton doublets of the SM by £ La and the RH charged leptons by ej? a , where a = 1, 2, 3 is the 
generation index. Moreover, we extend the SM by three RH neutrinos VR a which are total 
gauge singlets of the SM gauge group. 

Assume now that the early universe is baryon and lepton symmetric, i.e. the total 
baryon number B as well as the total lepton number L both vanish. Furthermore, there 
shall be no asymmetry between the LH and RH sectors. In other words, the baryon number 
B = Bl+Br and the lepton number L = Ll+Lr are both zero and vanish also in the LH and 
RH sectors, separately. Consider next the case where some process has produced a relative 
asymmetry L R = —L L > in the neutrino sector, i.e. L L and L R come from an excess of ii, a 
and URa, respectively. The LH and RH baryon and lepton numbers will be affected by two 
types of processes: (1) sphaleronic vacuum to vacuum transitions [14] and (2) left-right (LR) 
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Figure 2: LR equilibration processes £l + vr —> H and Ur £l + H changing Ll and Lr. 



equilibration processes (see Fig.|l|. Note that SU(2) L sphaleronic processes only act on the 
LH sector. They violate Bl and Ll by 3 units each, i.e. they are conserving Bl — Ll but not 
B L + L L - LR equilibration processes, by contrast, conserve B and L separately but violate 
B l ,r and Ll,r- Suppressing generation indices, the Yukawa couplings between 1 L and v R 
to the SM Higgs field H give rise to LR equilibration processes such as II + vr — > H and 
vr — > £l + H (see Fig. [2j. For sufficiently small Dirac neutrino Yukawa couplings, sphaleron 
processes will dominate the LR equilibration processes in the neutrino sector and Ll will be 
partly converted into a nonzero LH baryon asymmetry ABl > 0. However, LR equilibration 
processes in the baryonic sector, which also include SU(3) C sphaleronic processes [24], are 
fast and transfer half of ABl into a RH baryon number AB R > 0. Consequently, the 
SU(2)l sphaleronic processes reach equilibrium at 2Bl + Ll = 0. This changes Ll by an 
amount ALl = 2 ABl to a new lepton asymmetry L' L = Ll + ALl- Since Lr remains 
unaffected by sphalerons (L' R = Lr), the subsequent LR equilibration processes convert 
L' L and Lr into the final lepton asymmetries L" L and L" R , which are equal and given by 
L" LR = ALl/2 > 0. Thus, we arrive at a final total positive lepton asymmetry L" that 
is L" = L" L + L" R = ALl > 0. At the same time, the final total baryon asymmetry is 
B" = 2ABl = AB. As already mentioned, this requires the neutrino Yukawa couplings to 
be sufficiently small so as to allow the LR equilibration to take place only after sphaleron 
processes have dropped out of thermal equilibrium. A numerical study shows that this 
becomes possible for neutrino Yukawa couplings of the order < 10 -8 [9]. Therefore, if 
the neutrinos were Majorana particles, with the observed small neutrino mass scale m u ~ 
10 _1 eV generated by the type-I seesaw mechanism, the neutrino Yukawa couplings to the 
SM Higgs field would be ~ 1, which is by many orders too large. If, instead, the neutrinos 
are Dirac particles, the neutrino Yukawa couplings will be of the order m u /v ~ 1CT 12 , where 
v ~ 10 2 GeV is the SM Higgs vacuum expectation value (VEV), and sphaleron processes will 
dominate LR equilibration in the neutrino sector as required for successful Dirac leptogenesis. 

Let us briefly compare with the case of the SM (without massive neutrinos). In the SM, 
all Yukawa couplings in the quark and lepton sectors are ^> 10~ 8 and LR equilibration would 
be taking place roughly at the same time as the sphaleron processes. Starting with arbitrary 
B and L, the LR equilibration would then quickly drive B + L to values with B — L = 
such that sphaleron processes can only give B + L = and, thus, yield zero net baron and 
lepton asymmetries. To make the above mechanism for leptogenesis work, we thus require 
Dirac neutrinos which can provide sufficiently small neutrino Yukawa couplings. 
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3 Throat Geometry 



We will now be concerned with a model for Dirac neutrinos that generates the observed 
baryon asymmetry via Dirac leptogenesis as discussed in Sec. [2] Consider for this purpose 
three intervals in 5D flat space which are glued together at a single point as shown in Fig. [3] 
We will call the intervals throats. The coordinates on the three throats are respectively 
z\ l = yi), — (3^2/2), and z^ 1 = (x^,y 3 ), where the 5D Lorentz indices are denoted 
by capital Roman letters M — 0, 1, 2, 3, 5, while the usual 4D Lorentz indices are symbolized 
by Greek letters /i = 0,1,2,3. The coordinates yi,y 2 , and y 3 , describe the 5th dimension 
for the three throats. The physical space is thus defined by < y\ < nRi,0 < y 2 < vri? 2 , 
and < y 3 < vri? 3 , where Ri,R 2 , and R3, denote the size of the throats. The intersection 
point at yi = y% — y 3 = will be called the UV brane and the endpoints of the intervals 
at yi = irRi,y 2 = itR 2 , and y 3 = vri? 3 , will be denoted as IR branes. It will be useful to 
characterize this throat geometry by reflection symmetries interchanging the 1st and 2nd as 
well as the 2nd and 3rd throat (see Fig. [3]) . We will comment on these symmetries and how 
they are broken later. 

On the three throats, we assume the SM gauge group Gsm- The scalar sector contains 
three 5D Higgs doublet fields Hi,H 2 , and H 3 , that carry the same Gsm quantum numbers 
as the usual SM Higgs field. We assume that the three Higgs doublets live on separate 
throats: Each of the fields Hi,H 2 , and H 3 , is propagating on the first (Hi), second (-^2), 
and third (H 3 ) throat, respectively. By separating the scalar fields on the throats, we can, 
in the following, neglect the mixing between the scalars. The action of the Higgs doublets 
is then 

/ 3 pnRi 
/ dyi£ Hi , (1) 
i=i Jo 

where the 5D scalar Lagrangian density is 

C H . = (DmH^D" 1 ^ - HjH t + X l (HjH t ) 2 , (2) 

in which Dm is the covariant derivative. We will later discuss how the SM gauge group is 
spontaneously broken when only Hi acquires a non-zero VEV in the bulk of the 1st throat. 

Let us now focus on the lepton sector only (the discussion for quarks should be along 
the same lines). We suppose that the SM leptons propagate on all three throats. For this 
purpose, we start on each throat with 5D lepton fields L ita ,E ita , and N i>a , where i = 1,2,3 
labels the throat and a = 1, 2, 3 is the generation index. Each fermion with label i propagates, 
like Hi, only on the throat i. The fields L i a carry the quantum numbers of the LH SM lepton 
doublets, E i a of the RH SM charged leptons, while N i)CL are SM singlet neutrinos. As the 
5D action of the fermions we take 

/3 /* 7T R'l ' 

d 4 xJ2 / %(^ i ,ar M & i , a -m* a ^ !a ^ a ), (3) 

v=n,n, iv i,a=l ^° 

where \1/ = L, E, N, denotes the fermion species, mf a are the bulk masses of the 

Ki = Om - Ki, T M = (J 1 ® T 5 = diag(il 2 , -il 2 ), a = -1 2 , and j (i = 1, 2, 3) are the 
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Figure 3: Geometry of the three throats. The throats are described by intervals with coordi- 
nates yi, ?/2, and 2/3, that intersect at the UV brane at y\ — yi — y3 = 0. The Higgs doublets 
Hi, H 2) and H 3 , are living separately on the 1st, 2nd, and 3rd throat, respectively. All the 
other fields, such as the SM fermions, are propagating on all three throats. The throats are 
related by exchange symmetries D 1 and D 2 . 



Pauli matrices. Note that the fields L i a , E iA , and N ija , are vector-like in 5D. By imposing 
appropriate boundary conditions and interactions at the UV brane, we will later show how 
to obtain from these fields chiral fermion zero modes that can be identified with the SM 
fermions. 

We suppose that the 5D fermions couple to the scalar doublets only via Yukawa interac- 
tion terms localized at the IR branes of the throats: 

S Y = [ d 4 x V / Rl dy, t C 5 Yi + h.c, (4) 
J i=1 Jo 

in which the 5D Yukawa coupling Lagrangians C Yi are 

3 

4-* = 6 ^ - KRiWhabHih^ + {YXbH^E^} + h.c, (5) 

a, 6=1 

where Hi = io^i?*, while and are the complex 3x3 lepton Yukawa coupling matrices 
for Hi and Hi in the 5D theory. 
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We assume that the throats are subject to exchange symmetries (see Fig. [3]) which act 
on the scalar doublets Pj, the fermions ^fi >a , and the throat coordinates as 

D 1 : H x <-> H 2 , * lja «-> * 3ia , ^ <-> y 2 , (6) 

and 

D 2 : H 2 ^H 3 , tf 2)0 P£tf 3 ,6, 2/2 ^2/s- (7) 

Here, (Pj^) = P* are 3x3 matrix representations of some discrete flavor symmetry, i.e. the 
Pj, act on the generation indices. As a simple example, we will consider 




Is, (8) 



but other choices, such as those presented in the appendix, are also possible. The matrix 
P L yields a representation of an element of the group A (24) ~ (Z 2 x Z 2 ) x 5*3 (from the 



the symmetry D\ (D 2 ) 
ine between the throats 



class QC^) [25] and generates a Z 4 subgroup of A (24). In Fig. 
corresponds to a reflection symmetry with respect to the dashec. 
1 and 2 (2 and 3). Note that the symmetries D\ and D 2 require the throats to have equal 
lengths R\ = R 2 = P3. Moreover, for the choice of matrices in (J8|, the symmetries D\ and 
D 2 establish among the Yukawa coupling matrices the identities 

= Y 2 5 = p L r 3 5 , Yl = y 2 5 = P L Kf. (9) 

As we will see below, to obtain a realistic light fermion spectrum, we need to break D\ and 
D 2 at the UV brane. 

Let us briefly comment on how the symmetry D 2 could emerge from the product of two 
finite groups via spontaneous symmetry breaking. We begin with a group D' 2 ~ Z 2 xG, where 
G is the flavor symmetry group that acts on the fields on the third throat and is generated by 
the matrices in (|8]). The group Z 2 is given by the exchange symmetry Z 2 : X 2 <-> X 3 , y 2 «-> 
2/3, where X = H,L,E,N (generation indices have been neglected). Assume now a SM 
singlet scalar field S that carries a charge +1 under the Z 2 symmetry and transforms under 
application of the flavor symmetry transformation in ^ as S — > — S. When S acquires 
a nonzero VEV, D 2 will be broken at some high scale down to the subgroup D 2 of ([7]). 
However, instead of expanding further on the details of the possible underlying symmetry 
groups at high energies, we will, in the following, only be concerned with the symmetries Di 
and D 2 of the low-scale theory. 



4 Wavefunction Profiles 

In this section, we consider the boundary conditions for the 5D scalar and fermion fields, 
determine the mass spectra and wavefunctions of the bulk scalars, and demonstrate the 
exponential wavefunction localization of the fermion zero modes on the throats. 
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4.1 Scalar Boundary Conditions 

The scalar doublets on the three throats are supposed to be subject to the following BCs 



at the IR branes 
at the UV brane 



d yi H 1 



= 0, 
0, Ho 



3/1=0 "> - t - t 2,3 1 y 23 =o 



0. 



(10a) 
(10b) 



Note that we have on the first throat Neumann BCs at both endpoints for the field Hi, 
whereas the fields H 23 have Neumann BCs at the UV brane and Dirichlet BCs at the IR 
branes. The most general flat space KK expansion of the scalars, consistent with the BCs 



in (10a) is for Hi given by 
Hi(x^,yi) -- 



v^^Mcos^) 



n=l 



(11) 



while the KK expansions for the fields H 23 read 



Hi Ui 



(2n - l) Vi 




(12) 



n=l 



where i = 2,3. Note the important fact that the Dirichlet BCs at the UV brane have 
projected out the zero modes of H 23 , such that only Hi will have a zero mode. At the 
same time, the Neumann BCs at the IR branes ensure that the Hi are non-vanishing there. 
Di is broken by the different BCs for Hi and H 2 at the UV brane. Moreover, as we will 
explain further below, Di is broken by the bulk mass terms for the SM singlet neutrinos 
iVj i0 . Different from the symmetry Di, however, D 2 remains almost completely intact in the 
scalar sector. 

Denoting by M n (Hi) the mass of the nth KK state H\ n) of H at zero temperature, we 
thus arrive for Hi and H 2 3 at the mass squares of the KK states 



M 2 n (Hi) 



n 



7lR 1 



2-1 



M 2 n {H 



Hi 



2n - l\2i 



2nRi 



(i = 2,3), 



(13) 



where n — 0, 1, . . . for Hi and n — 1, 2, . . . for H 2 and if 3. Notice that the discrete symmetry 
D 2 in ((7]) establishes = jj-l as well as R 2 = R 3 . The mass squares M%(H 2 ) and M%(H 3 ) 
will therefore, up to small corrections, be practically degenerate. As we will see later, for 
successful Dirac leptogenesis, we will need small mass-squared splittings of the order 



\M 2 n (H 2 ) - M 2 n (H 3 )\/M 2 n (H 2 ) ~ 10- 



(14) 



Such small relative mass-squared splittings may be induced, e.g., at the quantum level, but 
we will not further specify an origin of this splitting here. The potential of the 5D scalar 



doublet Hi with the KK expansion given in (11) has a local minimum for the VEVs [26] 



wo 

V2 \v 



(Hi n) ) = 0, 



(15) 



where v is a real parameter with mass dimension +1 and n = 1,2,... In other words, 
only the zero mode acquires a non-zero VEV while all higher KK excitations have 
zero VEVs. Qualitatively, this is because only the zero mode has a negative mass square 
(coming from the potential), while the higher KK excitations have, for a sufficiently large 
compactification scale, always positive mass-squares. Similarly, since the zero modes of i?2,3 
have been projected out by the Dirichlet BCs, we will take for (n > 1) the VEVs 

= 0, (16) 
i.e. the VEVs of all KK excitations of H 2 ^ vanish. Therefore, only with the VEV given 



in (15) will be responsible for EWSB and for generating masses for the SM fermions from 



the Yukawa interactions in ([5]). 
4.2 Fermion Boundary Conditions 

Since fermions in 5D are vector-like, we have to impose appropriate BCs in order to obtain 
a chiral 4D theory. In doing so, we will apply the techniques introduced in [22,27] for multi- 
throat geometries. For this purpose, we write, neglecting generation indices, the 5D fermions 
on the ith throat (\1/ = L, E, N) as Dirac spinors of the form 

** = (t U ) > (17) 



Hi 



where ^ ta and denote two-component Weyl spinors and i = 1,2,3 labels the throat on 
which lives. The 5D action of the Dirac spinors is given by iS|, D in (J3|. In absence of 
brane-localized operators, the equations of motion for read 

-ia»d^ Li -d 5 V m + mfq m = 0, (18a) 
-i^d^ Ri + d^ Li + mf^ Ll = 0. (18b) 

Consider now for at both endpoints of the ith throat Dirichlet BCs ^Ri\o,nRi — 0, 
which lead for to the appearance of a single chiral zero mode with an exponential 5D 
wavefunction ~ exp(-mfyi) propagating on all three throats (for a discussion of exponential 
localizations of wavefunctions see [28,29]). This is achieved by connecting the throats by 
a brane-localized action 5uv at the UV brane. This action couples the to some extra 
fields ^ (i — 1,2) which are localized at the UV brane: 

Suv = [d A xJ2J2 r* dyi (Kv)^»?4 + h - c -) %). ( 19 ) 

J i=l j=l J ° 

where m^y is a UV brane mass parameter and /C* is a dimensionless 3x2 rank 2 matrix. 
Note that since the equations of motion in (18) are first order differential equations, the 



brane-localized operators will lead to a discontinuity of the wavefunction of the RH field 
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at yi = 0, i.e. lim e ^ ^Li(^, y% 



Ri(x tl ,y 



0) = 0. Including the brane interaction 



term (rrijjy) 2 fcfj ^{Ui) m (18a), we then obtain for the LH fields the BCs [22] 



E 

i=l 



/C.%a| w=o = 0. 



(20) 



As a consequence, two of the three zero modes decouple for large mjjy, leaving the remaining 
zero mode as a single chiral field propagating on all three throats. The wavefunction of this 
mode is 



A% exp(-mfyi)ipl(x) with Af 



r 3 

E 



1 — exp(— 2mf ttRi) 
2m? 



(21) 



where ip° L {x) is a 4D Weyl spinor with mass dimension 3/2. We can see from (21) that 
the actual localization of the zero modes depends on the signs of the bulk masses mf , 
which can be positive or negative. We thus obtain the following possible localizations: 
For mf 



rrin 



m{ 



mf < 0, the zero mode ^f L is localized at the three IR branes, for 
mf = mf > 0, it is localized at the UV brane, and for only one or two mf < 0, it is 
localized at the IR branes of the throats with positive mf . In a similar way, we can localize 
RH fermion zero modes on the UV and IR branes of the throats by replacing in the above 
considerations the LH and RH fields. 

The brane-localized Yukawa couplings at the IR branes in ^ lead to discontinuities of 
the wavefunctions of Liu, En, and Nu, at yi = nRi. The wavefunctions obey at the IR 
branes the BCs 





nR~ 




nR~ 




nR~ 


e for 


e — 



-OfHtNm + YfHiEK)^ 



tr; 



(22a) 
(22b) 
(22c) 



where = nRi — e for e — > (e > 0). Note that Lu,Eui, and Nm, are continuous 

over the whole interval, including both endpoints, i.e., in particular, E Ri \ 7rR - = E Ri \ vRi and 
^Ril-nRr = N Ri \ nR .. (In contrast to this, L Ri ,E Li , and N Li , are discontinuous at yi = nRi.) 

In our model, we have the LH and RH fermion zero modes L^\e^\ and with 
wavefunctions 



'R ' 



-.(0) 



yi 



E ( R '(x, yi ) 



A% exp(-mf yi ) £ L (x) 
Bq exp(mfyi)e R (x), 
exp(mfyi) u R (x), 



(23) 



that correspond, in the notation of Sec. |2j (up to a normalization) to £l and e R of the SM 
and to v R . Here, Bf denotes the normalization factor for the RH fields, which is similar to 
Af with only the signs in front of the bulk masses switched. We choose for the fermions 
the bulk masses as given in Tab. [T] As a consequence, the zero modes of the LH lepton 
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sgn(mf ) 


sgn(mf) 


sgn(mf) 


u 


-1 


-1 


-1 


Ei 


+1 


+ 1 


+ 1 




-1 


+ 1 


+ 1 



Table 1: Signs of the bulk masses for the fields Li,E{, and Ni, on the different throats. 




Figure 4: Wavefunction profiles of the scalars and fermions in the throat geometry 
(schematic). While the RH SM singlet neutrino zero modes iVj^ are localized towards 
the IR branes of the throats 2 and 3 (combined here as i — 2, 3), all SM fermion zero modes 
L/'jp and are symmetrically localized at all three IR branes of the throats. The Higgs 
zero mode lives only on the 1st throat giving a large overlap with all SM fermions but only 
a small overlap with the RH neutrino. The overlap of and the other fermions with the 
heavy scalars #23 j however, is large. 



doublets and the RH charged leptons E^ become symmetrically localized at the three 
IR branes of the throats. In contrast to this, the zero modes of the RH neutrinos are 
only localized towards the IR branes of the throats 2 and 3 but are repelled from the IR 
brane of the 1st throat. Schematically, the wavefunctions of the fields in the bulk, including 
the wavefunctions of the scalars, are depicted in Fig. [4 Note that and E^ have a 
large overlap with the zero mode which generates the lepton masses after acquiring a 
nonzero VEV. The overlap of the RH neutrinos with H[ ^ is, however, exponentially 
small, thereby suppressing the Dirac neutrino masses of the active neutrinos. The overlap 
of and with the higher KK-excitations H% and H^ n \ on the other hand, is larger, 
giving larger Dirac neutrino Yukawa couplings to the heavy scalars. We suppose for each 
particle species \I> that the bulk masses mf are flavor diagonal and degenerate for a fixed 
throat number i. This could, e.g., be ensured by an SU(3) flavor symmetry that is preserved 
on each throat but broken at the UV and IR branes. We will, however, not discuss further the 
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bulk flavor symmetry and its breaking but assume from now on simply the flavor-diagonal 
structure of the bulk masses and their degeneracy on each throat. 

In what follows, we will, for simplicity, go to the limit of strongly localized fermion zero 
modes, i.e. mfirRi ^> 1. In the low-energy effective theory, the lepton Yukawa couplings of 
the fermion zero modes are then described by the Lagrangian 

3 

C oS = Yl {(^l)^f 0) 4a^+(?i)a^ 0) ^,ae^ 
a,b=l 

+ [OQabHP£ L ,aVR, b + (fyabH^tLtaeRj,] } + h.C. + . . . , (24) 



i=2,3 



where Y{ oc Yf and Y{ oc Yf are the Yukawa coupling matrices of the 4D theory obtained after 



integrating out the extra dimension. In (24), we have included only the lightest KK scalars 
and restored the generation indices. The 4D Yukawa coupling matrices satisfy relations 
similar to those in Q: 

Y 1 <xY 2 = P L Y 3 , Y 1= Y 2 = P L Y 3 . (25) 

Since the symmetry D\ remains unbroken at the IR branes, the Yukawa coupling matrices 
Y\ and Y 2 ^ are, to leading order, in the 4D theory related by an overall rescaling factor 

Y x = ie- M l raW l Y 2 = * e-^l^l pLy 3i (26) 
v 2 v 2 

where R — R\ — R 2 — R% and \m N \ = \m^\ = \m 2 \ = \m^\. This connects directly 
the low-energy Yukawa coupling matrices Yx, accessible to neutrino oscillation experiments, 
with the Yukawa coupling matrices Y 2 $ that describe the interactions of the SM leptons with 
the heavy scalars ■ Calling the rescaling factor F = 4= exp (—2ttR \m N \) } we will later 
choose F ~ 10~ 5 to obtain realistic Yukawa couplings that give the right size ~ 10 _1 eV for 
the observed neutrino masses, while enabling, at the same time, successful leptogenesis. 
In <Suv, we wm assume the matrices /C* to be 




(27) 



where S <C 1 is a small symmetry breaking parameter. From (20), we thus have at the UV 



brane ^Li\ yi =o — ^ L2\ y2 =o = ^ ' L^\ y3 =o (for ^ Li as an example). We see that each field ^ u 



makes up 1/3 of the zero mode wavefunction. The matrices in (27) break the symmetry D\ 
but D 2 is only slightly broken. 

At a finite temperature T, the heavy scalars H± receive additive thermal corrections to 
their mass-squares, which are, neglecting Yukawa interactions, to leading order given by [30] 

AM„ 2 T (^) = -^j— T 2 , (28) 
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Figure 5: Decay amplitudes generating the decay asymmetries (i = 2,3). 

where X\ is the quartic self-coupling of the nth KK Higgs field on the ith throat. The 
symmetry D 2 establishes X^ = A3 , such that the corresponding leading order thermal 
corrections to the mass-squared splittings M^(H 2 ) and M^(H^) are zero. The fermion mass 
matrix /C* in (27) breaks D\ and D 2 . But since the Higgs fields Hi do not couple to the 



symmetry breaking terms at the UV brane, 5uv will only produce an unobservable shift in 
the potential without changing the thermal masses of the scalars [30] . 



5 Leptogenesis 

5.1 Bounds on Yukawa Couplings 

Let us now see how in our model the observed baryon asymmetry is generated via Dirac 
leptogenesis [9] through the decay of the heavy scalar doublets (i = 2,3). The scalars 
decay via 

H$$^e L + VR,h + e R , (29) 
involving the Yukawa coupling matrices Y 23 (see Fig. [5]). The allowed parameter space of 
the KK masses M n (H 2 $) and the corresponding neutrino Yukawa coupling matrices Y 2 ^ 
is restricted by several bounds. First of all, according to Sakharov's third condition, the 
asymmetry generating processes have to be out of equilibrium, i.e. in our scenario, the KK 
scalars have to decay at temperatures T < M 1 (H i ). Dirac leptogenesis, on the other hand, 
requires that all relevant decay processes have to take place at energies above the scale of 
EWSB T c , i.e. as long as the sphaleronic processes are in thermal equilibrium [14]. We thus 
have a time limit A£d for the relevant decays of the lowest scalar KK excitations 

A (D = t{T c) - tWH,)) = 0.30 x *| Qj - ^) , (30) 

where the time t is given by t — | H" 1 with the Hubble parameter if = 1.66 x i/gijf^ an d 
the factor denoting the number of accessible relativistic degrees of freedom at temperature 
T. 

The time A£d has to be compared with the life time t v for decays into RH neutrinos: 

T* - r; 1 = ^ . (31) 

tr(y^) M x {Hi) 
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Since the decays take place at temperatures that are small compared to the scalar masses, 
the Lorentz gamma factor can be neglected. For the asymmetry generating processes to take 
place we obviously need At D > r v . 

Moreover, for temperatures T > T c , in order to avoid LR equilibration by scattering 
processes, it is necessary that H(T) > Ts(T), where Ts(T) = n eq (as\v\) denotes the rate of ur 
annihilating scattering processes, n eq the equilibrium number density of target particles, and 
(cts|u|) is the thermally averaged scattering cross section. For H±, the dominant contribution 
from scattering processes at T c can be approximated as Ts(T c ) ~ A 2 (Y"i) 2 T c , where A ~ 1 is 
a typical gauge or Yukawa coupling. Consequently, the Yukawa couplings roughly satisy 



^ Vis ~ 10 ~ 8 ' (32) 

where we have neglected the generation indices. Observe that this condition is fulfilled by 
Dirac neutrinos with Yukawa couplings of the order Y\ ~ 10~ 12 that give neutrino masses 
~ 10 _1 eV consistent with observation. For the heavy scalars H^l-, the dominant scattering 
process takes the form [31] 

Ts^^^YifTl/MtiHi). (33) 
Setting the charged lepton Yukawa couplings Yj ~ 1, the bound on the Yj becomes 



unm (34) 

In order to measure the effectiveness of decays at T ~ M\(Hi), we introduce the quantity 



T(Hi) (Yf M, 



K = ^ = ^ ^P-. (35) 



2H{Mx(Hi)) 1.66 x 32^^ M 1 (H, 

For K <C 1, we are in a regime of pure "drift and decay" [31], called the weak wash-out 
regime, i.e. inverse decays are strongly suppressed and cannot erase the produced asymmetry. 

An upper limit on the KK masses M n (Hi) is set by the graviton bound [32,33], i.e. by the 
requirement that late decays should not spoil the production of light elements during Big 
Bang nucleosynthesis (BBN) at Tbbn ~ 0.3 MeV. On dimensional grounds, the KK graviton 
(G) decay rates take in our model the form 

r G( n) ~ n D 2 , (36) 

where ud denotes the number of decay channels and M n (G) is the mass of the nth graviton 
excitation. Note that due to the wavefunction profiles in our model the KK number is not 



conserved. From (36) we see that the KK graviton life time strongly depends on the KK 
masses, i.e. the compactification scale. Thus, we find two different parameter ranges: For 
Mi(G) > 50TeV, all KK excitations of the gravitons decay before the era of BBN and thus 
no further restrictions arise [32] . However, for the energy range we are particularly interested 
in, i.e. M\(G) < 50TeV, KK graviton decays lead to a bound on the reheating temperature 
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Figure 6: Summary of bounds on the Dirac neutrino Yukawa couplings 1^,3 to the heavy 
Higgs doublets H 2 and H 3 . The lines labeled by D, K, and LR, denote the exclusion regions 
for late decays (D), the weak washout regime (K), and LR equilibration (LR). Above T c , 
sphaleron processes are in thermal equilibrium and Trh is restricted below G by KK graviton 
decays. The preferred region considered here is Y 23 ~ 1CT 7 and Mi{H 2 ^) < 50TeV. 



T RH [33]. In this case, T RH must be close to the compactification scale and, consequently, 
only a few KK modes can come on-shell. In Fig. [6j we have summarized the bounds on the 
Dirac neutrino Yukawa couplings to the heavy Higgs doublets. The lines labeled by D, K, 
and LR, denote the exclusion regions for late decays (D), weak washout regime (K), and 
LR equilibration (LR). The allowed parameter space lies above the lines T c , D, K, and LR. 
The preferred region is therefore around Y 2 ,3 ^ 10~ 7 and Mi(F 2|3 ) ~ 1 TeV ... 50 TeVf] The 
neutrino Yukawa couplings to the heavy Higgs fields are therefore small but still by a factor 
~ 10 5 larger than the neutrino Yukawa couplings Y\ ~ 10~ 12 to the SM Higgs h[°\ generating 
the observed neutrino masses. In the flat limit, the mass range for M\{H 2 ^) translates into 
the range R^ 1 ~ 1... 50 TeV for the compactification scale, implying a fundamental Planck 
scale of the order M* = ~ 10 13 . . . 10 14 GeV. 

5.2 Lepton Asymmetry 

The lepton asymmetry is generated by the interference between the tree-level and one- 
loop amplitudes shown in Fig. [5] The one-loop amplitude must involve on-shell LH lepton 
doublets and RH charged leptons. The decay of the KK mode , e.g., leads to the decay 

1 For Y 2 ~ 10 -12 of the order the Dirac Yukawa couplings of the active neutrinos, we would need to go to 
a resonant limit with very tiny relative scalar mass-squared splittings ~ 1CP 20 . 
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asymmetry [34,35] 

» 



V{H { 2 n) -> i x v R ) - T(Hj n) -> EjUr) 
lm[tr(Y}Y 3 )tr(Y]Y 3 )} M 2 n {H 2 



(37) 



87rtr(Y 2 t y 2 ) M%(H 3 ) - M*(H 2 )' 
where we have used the fact that the dominant contribution comes from the pair of KK states 



with the same level number n [36]. In (37), we have only taken the self-energy contributions 



into account, since additional vertex corrections can be neglected in the resonant limit. Note 



that our expression for the resonantly enhanced lepton asymmetry in (37) corrects the result 
in [9] in two ways: We have (i) included a complex conjugation of the Yukawa coupling 
matrix Y 2 and have (ii) taken in the numerator the product of two traces instead of a single 
trace. 



Equation (37) holds as long as M n (H 3 ) — M n (H 2 ) ^> T n (H 2 ), which is satisfied for our 
range of parameters with small Yukawa couplings. Similarly, the decay of H 3 leads to a 



decay asymmetry e 3 which is obtained from the expression for by interchanging in (37) 



the Yukawa coupling matrices Y 2 <-> Y 3 and the Higgs fields H 2 <-> H 3 . For (Y 2 ) a b, (Y 3 ) a b < 
1CT 7 , a scalar mass-squared splitting of the order ~ 1CT 8 as given in (14) produces a total 
decay asymmetry e^ al = + £3 of the order e^ al ~ — 1CT 8 . For this choice of parameters, 
an equal number n„(i) of and h!> (i = 2, 3) leads to a net number density in the RH 

i 

neutrino sector ni£ — 4 e^lai 71 „i» . The out-of equilibrium decays in the "drift and decay" 

i 

limit [31], i.e. n ( n ) ~ n 7 , thus yield a neutrino number to entropy ratio 

i 

n max (n) »m» , (n) 

F " = £ 5t ~ £ ^ < 38 > 

n=l n=l 9* 



In (38), the dominant contribution to the asymmetry is generated around the energy at 



which #23 drop out of thermal equilibrium such that the entropy is given by = g^n^, 

where gi n ^ is the number of relativistic degrees of freedom at T ~ M n (H 23 ). In the standard 
model, we have ~ 100. This is altered in the 5D model by taking the additional KK 
excitations into account, giving = g* + gf K Q(T — i?~ 1 )Ti?, where Q(x) = 1, for x > 0, 
and O(x) = 0, for x < 0. The number of relativistic degrees of freedom for the zero modes 
is gi°^ ~ 100 and we take also gf K ~ gi°\ The asymmetry e(™ tal is independent of n, since 
the scalar mass-squared splittings as well as the expression for the Yukawa couplings are 
independent from n. Thus, we can approximate the final asymmetry by 



'tot a 

g* 



Y v ~ ^tai _ _ 10 -i°. (39) 



2 For a large number n max of scalar KK excitations the sum over 1 /gi can be approximated logarithmi- 
cally [36] by J2n=T(9*°^ + n 5? K ) _1 ~ m ( g *gJ+^^* ) ano " wc * nus obtain a factor 0(1) for a large range 
of energies instead of a factor ~ 1CP 2 , as in the case here of only a few KK states. 
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An analysis of chemical potentials [37] reveals that for initial B — L = and all the heavy 
KK excitations decaying out of equilibrium, the number densities of baryons and leptons are 
related to the number density of RH neutrinos by n B = Ul = — f§ n UR [9] . This means that 
Y v is converted by sphaleron processes into a baryon asymmetry Y B ~ —Y v of the order the 
observed value Y B = (8.62 ± 0.27) x KT 11 [38]. 



5.3 Correlation of Low-Energy Parameters 

The Yukawa coupling matrices appearing in the decay asymmetries e% and exhibit the 
important feature that they are related to the Yukawa couplings of the low-energy theory 
by the exchange symmetries D\ and D2 in ^ and ([7]). The BAU in our model becomes 
therefore connected with the low-energy neutrino masses, mixing angles, and the Dirac 
CP phase observable in neutrino oscillations. In the basis where the charged lepton mass 



matrix is diagonal, the trace over the neutrino Yukawa couplings in (37) takes the form 
tr(Y 2 ^F 3 ) = tr[L^PMNS^ diag ^3 diaS ^pMNS-^ >it ]- Inserting the Yukawa couplings in (25) into the 
expression for e^l, one can then study the BAU as a function of the solar, atmospheric, and 
reactor mixing angles #12, #23, and #13, and the Dirac CP phase 5, of the low-energy leptonic 
mixing matrix ?7pmns r mi ' 



In Fig. [7| we show the correlations between the reactor angle #13, the atmospheric angle 
823, and the Dirac CP phase 5, for the matrix representation P L in ([8j that generates a Z A 
subgroup of A(24). The correlations are given as a function of the Higgs mass (-squared) 



splitting in (14) for the lowest KK excitation (n = 1). The solar angle and BAU are fixed 
at their best-fit values # 12 = 33.2° and Y B = 8.62 x 10~ n . We have assumed a normal 
hierarchical neutrino mass spectrum with mass ratios mi : '■ m 3 = 0.04 : 0.2 : 1 (similar 
results can, however, also be obtained for an inverted neutrino mass hierarchy). Denoting by 
Yj diag the matrix obtained after diagonalization of the 4D Dirac neutrino Yukawa coupling 
matrix Yi, we have, in Fig. [7, taken (Y 2 diag )33 = (Y 3 diag )33 = 10~ 7 , while the mass splitting 
of the decaying scalars has been varied in the range 0...10 -8 . The region for 5 > ir is 
unphysical, since it would lead to a wrong sing of the BAU. In the appendix, we present 
the results for the same parameters with P L taken as a matrix representation of elements 
of the groups A (27) and A (54). The correlations between the low-energy lepton mass and 
mixing parameters make our model testable at future neutrino oscillation experiments such 
as Double Chooz [41], T2HK [42], or a neutrino factory [43]. 



6 Summary and Conclusions 

In this paper, we have presented a model for resonant Dirac leptogenesis on a 5D flat multi- 
throat background. The baryon asymmetry is generated by the decay of heavy scalars, which 
are copies of the SM Higgs. The throats which are subject to discrete exchange symmetries 
allow to solve several possible shortcomings of the original scenario for Dirac leptogenesis. 
First, the model provides an origin of the heavy decaying scalars as KK excitations of 5D 

3 For a discussion of possible forms of Dirac neutrino mass matrices see also [40] . 
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Figure 7: Correlations between 6*13 and 5 (a), 623 and 5 (b), and #13 and #23 (c), for P L in 
(Jsj) as a function of the relative Higgs mass(-squared) splitting. In (a), (b), and (c), we have 
assumed the values 623 = 45.0°, 6> 13 = 9.1°, and 5 = n/2, respectively. 



Higgs fields. Second, the exchange symmetries protect a near mass degeneracy of the scalars 
which leads to resonant decays. This enables Dirac leptogenesis at energy scales as low as 
1 TeV ... 50 TeV that may be in reach of a collider. Third, the discrete symmetries, which 
are broken in the bulk, connect the observed BAU with the Yukawa couplings of the low- 
energy theory. This leads in our model to non-trivial correlations between the lepton mixing 
parameters. We have studied the dependence of the BAU on the atmospheric angle, the 
reactor angle, and the Dirac CP phase for several discrete group representations and found 
strong correlations between the mixing angles and the CP phase. This makes our model 
testable at future neutrino oscillation experiments such as neutrino factories. 

It would be interesting, e.g., to consider in more detail the Boltzmann equations for our 
model, to study the throat stabilization necessary to ensure the resonant decays, and to 
investigate possible collider implications. 
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A Correlations for Elements of A (27) and A (54) 



Let us now present further results for the correlations between the low-energy lepton mixing 
parameters, when the matrix P L in ^ is a matrix representation of a generator of a cyclic 
subgroup of the groups A (27) or A (54) (for a discussion of A (54) as a flavor symmetry, 
see [44]). In what follows, we will, as in Fig. [7j set throughout the solar angle and the BAU 
equal to their best fit values 6 12 = 33.2° and Yb = 8.62 x 10~ n . Moreover, the neutrino 
Yukawa couplings are chosen as in Sec. 5^3 As a first example, let us consider for P L the 
following matrix representation: 




(40) 



where u = exp (2vri/3). The matrix P L is taken from the class 3 C[ 1,1] [45] or C 9 [46] of A (27) 
and generates a Z3 symmetry. Fig. [8] shows the correlation between #13 and 5 for this choice 
of P L . The Higgs mass (-squared) splitting is defined as for Fig. uj Note in Fig. Mthat we 
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Figure 8: Correlation between the reactor angle #13 and the CP phase S, for 623 = 42.1° and 



the matrix P L in (40) 



have taken for the atmospheric angle the value 623 = 42.1° (for 823 = 45° the BAU would 
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become further suppressed by two orders of magnitude due to an accidental cancelation). 
We show in Fig. [8] only the correlation between # 13 and 8 since it is much stronger than the 
dependencies of these parameters on 6 I 23 . Radically different results are obtained for another 
group element from the same class of A (27) with matrix representation 

(uo 2 0\ 

P L = uj . (41) 
\0 l) 

This matrix leads to the correlations between low-energy parameters shown in Fig. |9j We 
observe that in this case there is no strong dependence of # 13 on 8. 
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Figure 9: Correlations between #13 and 8 (a) and between #13 and #23 (b) for P L as in (41) 
In (a) and (b) we have set 823 = 45.0° and 8 = n/2, respectively. 



As a final example, consider the matrix representation 

/0 w 0\ 

P L = 1 (42) 
\0 to 2 J 

for an element taken from the class 9C^ [25] of A (54), which generates a Zq symmetry. 



In Fig. (10), we have summarized the resulting correlations between # 13 , 9 2 s, and 8, for this 
P L . The Higgs mass(-squared) splitting is defined as for the other examples. We see that 
the correlations between the leptonic mixing parameters differ strongly from those shown in 
the other figures. 
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Figure 10: Correlations between 9 13 and 5 (a), #23 and 5 (b), and 9\ 3 and #23 (c) 
matrix P L in Q. In (a),(b), and (c), we have set 9 23 = 45.0°, 13 = 9.1°, and 5 
respectively. 
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